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Abstract 

Using data from the FOCUS (E831) experiment, we have searched for T violation 
in charm meson decays using the four-body decay channels D° — > K~ K + ir~7r + , 
D + — > KgK + ir~ir + , and Df — > KgK + ir~7r + . The T violation asymmetry is ob- 
tained using triple-product correlations and assuming the validity of the CPT the- 
orem. We find the asymmetry values to be A^ vio i(D°) = 0.010 ± 0.057(stat.) ± 
0.037(syst.), A Tvio i(D + ) = 0.023 ± 0.062(stat.) ± 0.022(syst.), and A Tvioi (D+) = 
—0.036 ± 0.067(stat.) ± 0.023(syst.). Each measurement is consistent with no T vi- 
olation. New measurements of the CP asymmetries for some of these decay modes 
are also presented. 



1. Introduction 

The origin of CP violation remains one of the most important open questions 
in particle physics. Within the Standard Model, CP violation arises due to 
the presence of a phase in the Cabibbo-Kobayashi-Maskawa (CKM) quark 
mixing matrix. Although the main focus has been on rate asymmetries, there is 
another type of CP violating signal which could potentially reveal the presence 
of physics beyond the Standard Model. Triple-product correlations of the form 

1 See http://www-focus.fnal.gov/authors.html for additional author information. 
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v[ ■ (vq, X V3), where each vl is a spin or momentum, are odd under time reversal 
(T). By the CPT theorem, a nonzero value for these correlations would also 
be a signal of CP violation. A nonzero triple-product correlation is evidenced 
by a nonzero value of the asymmetry [1] 

= r>i • (v 2 x v 3 ) > 0) - r>i • (v 2 x v 3 ) < 0) 
T ~ r(«i ■ (tr 2 x tr 3 ) > o) + r(vi ■ (? 2 x «i) < o) u 

where T is the decay rate for the process. There is a well-known technical 
complication: strong phases can produce a nonzero value of At, even if the 
weak phases are zero, that is CP and T violation are not necessarily present. 
Thus, strictly speaking, the asymmetry At is not in fact a T- violating effect. 
Nevertheless, one can still obtain a true T-violating signal by measuring a 
nonzero value of 



^Tvioi = \{A T - A T ) (2) 

where At is the T-odd asymmetry measured in the CP-conjugate decay pro- 
cess [2]. 

This study was inspired by a paper of Ikaros Bigi [3]. In this paper Bigi 
suggested a search for T violation by looking at the triple-product corre- 
lation (using the momenta of the final state particles) in the decay mode 
D° —> K~ K + 7r~7r + . Such a correlation must necessarily involve at least four 
final-state particles. This can be understood by considering the rest frame of 
the decaying particle and invoking momentum conservation. The number of 
independent three-momenta is one less than the number of final-state parti- 
cles, so a triple product composed entirely of momenta requires four particles 
in the final state [4]. 

We calculate A Tviol for the decay modes D° — > K~ K + n~iv + and Df, — > 
KgK + 7T~7T + using data from the FOCUS experiment. 



FOCUS is a charm photoproduction experiment [5] which collected data dur- 
ing the 1996-97 fixed target run at Fermilab. Electron and positron beams 
(with typically 300 GeV endpoint energy) obtained from the 800 GeV Teva- 
tron proton beam produce, by means of bremsstrahlung, a photon beam which 
interacts with a segmented BeO target. The mean photon energy for triggered 
events is ~ 180 GeV. A system of three multicell threshold Cerenkov coun- 
ters performs the charged particle identification, separating kaons from pions 
up to 60 GeV/c of momentum. Two systems of silicon microvertex detectors 
are used to track particles: the first system consists of 4 planes of microstrips 
interleaved with the experimental target [6] and the second system consists 
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of 12 planes of microstrips located downstream of the target. These detec- 
tors provide high resolution in the transverse plane (approximately 9 /xm), 
allowing the identification and separation of the primary (production) and 
the charm secondary (decay) vertices. Charged particle momentum is deter- 
mined by measuring deflections in two magnets of opposite polarity through 
five stations of multiwire proportional chambers. 

2. Search for T violation in the decay mode D° — > K~K + ir~ir + 

The decay mode D° — > K~ K + 7T~ir + is Cabibbo-suppressed and may be pro- 
duced as a non-resonant final state or via two-body and three-body intermedi- 
ate resonant states. In a previous paper we determined its resonant substruc- 
ture and the branching ratio T(D° -> K~ K + n-n + )/Y(D Q -> K-ti-tx + ti + ) [7]. 

The final states are selected using a candidate driven vertex algorithm [5]. A 
secondary vertex is formed from the four candidate tracks. The momentum 
vector of the resultant D° candidate is used as a seed track to intersect the 
other reconstructed tracks and to search for a primary vertex. The confidence 
levels of both vertices are required to be greater than 1%. Once the production 
and decay vertices are determined, the distance L between the vertices and 
its error cxl are computed. The quantity L / <7l is an unbiased measure of the 
significance of detachment between the primary and secondary vertices. This 
is the most important variable for separating charm events from non-charm 
prompt backgrounds. Signal quality is further enhanced by cutting on Iso2, 
which is the confidence level that other tracks in the event might be associated 
with the secondary vertex. We use L / o L > 6 and Iso2 < 10%. We also require 
the D° momentum to be in the range 25-250 GeV/c (a very loose cut) and the 
primary vertex to be formed with at least two reconstructed tracks in addition 
to the D° seed. 

The Cerenkov identification cuts used in FOCUS are based on likelihood ratios 
between the various particle identification hypotheses. These likelihoods are 
computed for a given track from the observed firing response (on or off) of all 
the cells that are within the track's {(3 = 1) Cerenkov cone for each of our three 
Cerenkov counters. The product of all firing probabilities for all the cells within 
the three Cerenkov cones produces a x 2 -like variable Wi = —2 ln(Likelihood) 
where i ranges over the electron, pion, kaon and proton hypotheses [8]. All 
kaon tracks are required to have Ak = W n — Wk (kaonicity) greater than 3; 
whereas all the pion tracks are required to be separated by less than 5 units 
from the best hypothesis, that is picon = W min — W n (pion consistency) is 
greater than —5. 

In addition to these cuts (also used in our previous analysis of this decay 
mode), we require a _D*-tag. The sign of the bachelor pion in the D*^ decay 
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chain D*+H 7r +< - ) is used to identify the neutral D as either a D° 

or a £)°. We require that the mass difference between the D° and the D* mass 
be within 4 MeV/c 2 of the nominal mass difference [9]. 

Using the set of selection cuts just described, we obtain the invariant mass 
distributions for K~K + tt~tt + shown in Fig. 1, where the first plot is the total 
sample and the other two plots show the D° and D° samples separately. 

The mass plots are fit with a function that includes two Gaussians with the 
same mean but different sigmas to take into account different momentum res- 
olutions in our spectrometer [5] and a second-order polynomial for the combi- 
natorial background. A log-likelihood fit gives a signal of 828±46 K~ K +/ k~tt + 
events for the total sample, 362 ± 31 D° events, and 472 ± 34 D° events. The 
fitted D° masses are in good agreement with the world average [9] and the 
widths are in good agreement with those of our Monte Carlo simulation. 

From the D° sample we can form a T-odd correlation with the momenta: 

C T = PK+ ■ (Pn+ X Pit-) (3) 

and from the D° sample we form: 

C T = PK- ■ (Ptt- *Pit+)- (4) 

As we have seen in the introduction, finding a distribution of Ct different from 
— Ct establishes CP violation [3]. 

Fig. 2 shows D°(D°) signals separated by the sign of Ct{Ct)- A log-likelihood 
fit , with the same fit function described previously, gives the yields summarized 
in Table 1. 
Table^ 

D° (D°) yields split by C T (C^) sign. 



Decay mode 


Request 


Events 


D° -> K~K + 1T~TT + 

D° -» K'K+tt-tt + 

~W -> K~K + 1T~1T + 


C T > 
C T < 
CV > 
CV < 


174 ± 21 
190 ± 24 
255 ± 24 
220 ± 25 



Before forming the asymmetry At (At) we have to correct for detection ef- 
ficiencies, accounting for possible differences in spectrometer acceptance and 
Cerenkov identification efficiency for positive/negative kaons and pions. 2 This 

2 It is well-known that in fixed-target experiments there are production asymme- 



5 



is, however, a small effect. From the efficiency corrected yields we compute the 
asymmetry: 



_ r(c T > o) - r(c T < o) 
T t(c t > o) + r(c T < o) 



(5) 



and 



i^-^>o)-r(-^<o) 

T(-C T > 0) + T(-C T < 0)' V ; 



The resulting T-violation asymmetry Axvioi is: 



A Tviol = -(A T - A T ) = 0.010 ±0.057. (7) 

Without the efficiency correction it would have been A^moi = 0.014 ± 0.057. 

This determination has been tested by modifying each of the vertex and 
Cerenkov cuts individually. Although the statistics is limited, the T-violation 
asymmetry is stable versus several sets of cuts as shown in Fig. 3. All the 
measurements are consistent with for the T-violation asymmetry. 



3. Search for T violation in the decay mode D — > KgK + ir n + 

The decay channel D + — > K s K + tt~tt + is Cabibbo-suppressed and like D° — > 
K~ K + 7T~ir + , it may be produced as a non-resonant final state or via two- 
body and three-body intermediate resonant states. Its relative branching ratio 
T(D + -> K° s K + Tr~TT + )/T(D + -> K° s TT-n + 7r + ) has been measured [10]. TJ+ -> 
KgK + 7T~7T + is observed in the same histogram as D + — > KgK + ir~iT + and we 
fit for both signals. 

The final states are selected using a candidate driven vertex algorithm as de- 
scribed in the previous section. The Kg is reconstructed using techniques 
described elsewhere [11]. The K s and the charged tracks are used to form a 
D candidate which is used as a seed track to intersect the other reconstructed 
tracks and to search for a primary vertex. The confidence levels of both ver- 
tices must be greater than 1%. We also use L / a L > 6 and Iso2 < 1% and 
require the primary vertex to be composed of at least two reconstructed tracks 
in addition to the D seed. 

tries between charm and anticharm particles. As a result the D° momentum distri- 
bution is different from the D° distribution. 
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Using these selection cuts, we obtain the invariant mass distributions for 
KgK + 7T~7T + shown in Fig. 4, where the top plot is the total sample and the 
bottom two plots show the D and D samples separately. 

The mass plots are fit with a function that includes a Gaussian for the D + and 
a Gaussian for the with the widths fixed to those given from our Monte 
Carlo simulations. We use a second-order polynomial for the combinatorial 
background in addition to two reflection peaks from A+ — > pKgir~7T + and 
D + — > A'^7r + 7r~7r + . The A+ yield is fixed after first fitting the sample with 
the K + mass changed to the proton mass. The D + — > A^7r + 7r _ 7r + yield is 
determined by using the Monte Carlo misidentification rate of a pion as a 
kaon and the yield of D + — > KgK + 7r _ 7r + . A log-likelihood fit gives a signal 
of 523 ± 32 events for the and a signal of 508 ± 34 events for Df. The 
K° s K + ti-ti + sample has 240 ± 22 D + and 270 ± 25 D+ events, while the 
K° s K-^-^ + sample has 282 ± 23 D~ and 239 ± 24 Dj events. The fitted 
D masses are in good agreement with the world average [9]. Also the excess 
of D~ over D + events is consistent with more D° mesons than D° mesons 
being produced. These photoproduced excesses have been observed in previous 
higher statistics studies by FOCUS [12,13]. 

The mass plots shown in Fig. 5 are the Df-> (D7 S A signals split by the sign of 
Ct{Ct)- A log-likelihood fit, with the same fit function described previously, 
gives the yields summarized in Table 2. 

Table 2 

Dt s) (D7 s) ) yields split by C T (C T ) sign. 



Final State 


Request 


D + Events 


Df Events 


K%K+ir~ir + 
K° s K--k-tt + 


C T > 
C T < 
Ct~> 
Ct~< 


122 ± 16 
118 ± 16 
145 ± 16 
137 ± 16 


126 ± 17 
147 ± 18 
120 ±17 
119 ±16 



After correcting for detection and reconstruction efficiencies as given by the 
Monte Carlo simulation, we form the asymmetry At (At) as given by equation 
(5) and equation (6): 

A Tvi oi(D + ) = ^(A T - At) = 0.023 ± 0.062 (8) 

AiwW) = \(A T -1W) = -0.036 ± 0.067. (9) 
Without the efficiency corrections the numbers are essentially the same. A 
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scan of ^Tvioi under a variety of different selection criteria is presented in 
Fig. 6. 

4. Systematic Uncertainties 

Systematic uncertainties on the T- violation asymmetry measurement can come 
from different sources. We determine five independent contributions to the 
systematic uncertainty: the split sample component, the fit variant compo- 
nent, the component due to the particular choice of the vertex and Cerenkov 
cuts (discussed previously), the dilution due to an erroneous D* tag for the 
D° — > K~ K + 7t~tt + channel, and a component due to the limited statistics of 
the Monte Carlo. 

The split sample component addresses the systematics introduced by a resid- 
ual difference between data and Monte Carlo, due to a possible mismatch in 
the reproduction of the D momentum and the changing experimental condi- 
tions during data collection. This component has been determined by splitting 
data into four independent subsamples, according to the D momentum range 
(high and low momentum) and the configuration of the vertex detector, that 
is, before and after the insertion of an upstream silicon system. A technique, 
employed in FOCUS and in the predecessor experiment E687, modeled after 
the S-f actor method from the Particle Data Group [9] , is used to try to separate 
true systematic variations from statistical fluctuations. The T-violation asym- 
metry is evaluated for each of the 4 (= 2 2 ) statistically independent subsamples 
and a scaled variance a (that is the errors are boosted when \ 2 / (N — 1) > 1) is 
calculated. The split sample variance cr split is defined as the difference between 
the reported statistical variance and the scaled variance, if the scaled variance 
exceeds the statistical variance [14]. 

Another possible source of systematic uncertainty is the fit variant. This com- 
ponent is computed by varying, in a reasonable manner, the fitting conditions 
on the whole data set. In our study of the D° mode, we fixed the widths 
of the Gaussians to the values obtained by the Monte Carlo simulation, we 
changed the background parametrization (varying the degree of the polyno- 
mial), we modified the fit function in order to take into account the reflection 
peak from D° — » K~7t + 7t~tt + [7], and we use one Gaussian instead of two. 
For all modes, the variation of the computed efficiencies due to the differ- 
ent resonant substructure simulated in the Monte Carlo has been taken into 
account. The T-violation values obtained by these variants are all a priori 
equally likely, therefore this uncertainty can be estimated by the r.m.s. of the 
measurements [14]. 

Analogously to the fit variant, the cut component is estimated using the stan- 
dard deviation of the several sets of cuts shown in Fig. 3 and Fig. 6. Actually, 
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this is an overestimate of the cut component because the statistics of the cut 
samples are different. 

An erroneous D* tag can obviously dilute the measured asymmetry Aivioi- We 
find a dilution 3 of 0.9846 ± 0.0029 for the D° sample and 0.9882 ± 0.0025 for 
the D° events (the dilutions are slightly different because, as we have already 
seen, the D° and D° momentum distributions are different). Then we com- 
puted the ^Tvioi asymmetry taking into account this dilution and estimated 
the uncertainty by using the difference between this determination and the 
standard one. 

Finally, there is a further contribution due to the limited statistics of the Monte 
Carlo simulation used to determine the efficiencies. Adding in quadrature all of 
these components, we obtain the final systematic errors which are summarized 
in Table 3. 

Table 3 



Contribution to the systematic uncertainties of the T- violation parameters for D , 
D+, and D+ 



Source 


D° Uncertainty 


D + Uncertainty 


Df Uncertainty 


Split sample 


0.000 


0.000 


0.000 


Fit Variant 


0.009 


0.006 


0.004 


Set of cuts 


0.035 


0.021 


0.022 


.D*-tag dilution 


0.002 






MC statistics 


0.009 


0.004 


0.006 


Total systematic error 


0.037 


0.022 


0.023 



5. Conclusions 

Using data from the FOCUS (E831) experiment at Fermilab we have searched 
for T violation in charm meson decays. It is a clean and alternative way to 
search for CP violation. This is the first time such a measurement has been 
performed in the charm sector. 

We determine the final values for the T-violation asymmetries to be 
^TvioiP ) = 0.010 ± 0.057(stat.) ± 0.037(syst.), 



3 The dilution, measured by means of our Monte Carlo simulation, is defined as 
D = ^T_jy ■ For the D° (D°) dilution, R is the number of generated D° (D°) events 
reconstructed correctly as D° (D°) and W is the number of generated D° (D°) 
events reconstructed as D° (D°). 
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A Tviol (D+) = 0.023 ± 0.062(stat.) ± 0.022(syst.), and 
A Tviol (D+) = -0.036 ± 0.067(stat.) ± 0.023(syst.). 

It is interesting to compare the Axvioi measurements with the usual CP asym- 
metry measurements. Following the procedure described in a previous pa- 
per [12] we determine Aqp for D° —>■ K~ K + ir~ir + and D + — > KgK + 7i^n + , 
where we used D° K~7r + 7T^7i + and D + — > Kg7i + 7i^7i + to account for dif- 
ferences at production level. Systematic errors are obtained from the same 
sources and in the same manner as for the At V wI measurement. 

We measure: 

A CP (D°) = -0.082 ± 0.056(stat.) ± 0.047(syst.) and 
A CP (D + ) = -0.042 ± 0.064(stat.) ± 0.022(syst.). 

Both ^Tvioi and Aqp are consistent with zero. While our measurements are 
consistent with no T violation, we encourage higher statistics experiments to 
repeat these measurements. 
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Fig. 1. K K + tt 7r + invariant mass distributions for: (a) total sample 
D*+(-) D°(I)0)7r + (-), (b) £>° sample, £>*+ -> £> tt + and (c) W sample, 
D*~ — > D°tt~ . The fit (solid curve) is explained in the text. 
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Fig. 3. T- violating asymmetry Ax v ioi versus several sets of cuts. We varied the 
confidence level of the secondary vertex from 1% to 5% (5 points), Iso2 from 10~ 6 
to 1 (7 points), L / ctl from 5 to 15 (11 points), Ak from 1 to 5 (9 points), picon 
from —6 to —2 (9 points). The dashed lines show the quoted At viol asymmetry ±1<t. 
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Fig. 5. Upper plots are invariant mass plots for D + — > KgK + Tr + -K~ for Ct > 
in Fig. 5a and for Ct < in Fig. 5b. Lower plots are invariant mass plots for 
D~ — > K^K~it + tt~ for Ct > in Fig. 5c and for Ct < in Fig. 5d. All plots have 
an L/a > 6 cut. The shaded regions shown in each plot are explained in the text 
and in Fig. 4. 
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Fig. 6. T-violating asymmetry ^Tvioi versus several sets of cuts. We varied the ul 
from 4 to 12 (8 points), Iso2 from 10~ 2 to 10~ 6 (5 points), confidence level of the 
secondary vertex from 1% to 5% (5 points), from 2 to 5 (7 points). The dashed 
lines show the quoted ^Tvioi asymmetry ztlu. 
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